
function dx=fun(t,x,u)
global Q1 Q2 Q3
Q1=u(1);
Q2=u(2);
Q3=u(3);

f1=Q1*x(1)+x(2)*x(3);
f2=-10*x(2)+Q2*x(3);
f3=-x(1)*x(2)+Q3*x(2)-x(3);
dx=[f1;f2;f3];
end